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METHODS AND SYSTEMS FOR ANALYZING FLUTTER TEST DATA USING 
10 NON-LINEAR TRANSFER FUNCTION FREQUENCY RESPONSE FITTING 

Inventor 
Charles E. Goodman 

Field of the Invention 
15 The present invention relates generally to analyzing data and, more specifically, to 

analyzing flutter test data. 



Background of the Invention 
In the field of aeroelasticity, the phenomenon known as flutter generally refers to a 
20 condition produced by the coalescing and proper phasing of two or more structural vibration 
modes of a structure, such as an aircraft wing, fuselage, empennage, or other stractural 
component. During flight, the aerodynamic forces over such structures may cause an 
unstable oscillatory aeroelastic deformation of the structure referred to as flutter. Flutter of 
an aeroelastic structure may depend on numerous factors, including the mass, stiffiiess, and 
25 shape of the structure, and the particular operating conditions of the structure, including the 
velocity and density of the airstream. 

Flutter of an aircraft wing may involve both bending and torsional types of motion. 
In some cases, the aeroelastic deformations may be relatively mild and stable within the 
normal operating envelope of the aircraft. In the case of flutter, however, the aeroelastic 
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deformations are driven into an unstable mode in which the torsional motion extracts energy 
from the airstream and drives the bending mode to increasingly higher amplitudes, causing 
oscillations of increasing amplitude that eventually result in catastrophic failure of the 
structure. 



maximum allowable flight speed before flutter is encountered are critical priorities for 
designers of aeroelastic structures and aerospace vehicles. Exhaustive flight and wind timnel 
tests are usually conducted to record and observe the flutter characteristics of the various 
aeroelastic structures of an aircraft over the entire operating envelope of the aircraft. For 

10 military aircraft, the attachment of external stores to the aircraft fiirther complicates the 
analysis and increases the extent and complexity of the testing due to the multiple store 
configurations dictated by modem combat, surveillance, and fiiel capacity requirements. 
These differences in store number, type, and mounting location give rise to complex multi- 
variable oscillation coupling pattems, and can give an aircraft as many different flutter 

15 speeds as there are store configurations. 

The acquisition and analysis of flutter test data may be difficult and expensive, 
particularly when considerable flight testing is involved. The two principal methods of 
acquiring flutter flight test data are so-called dwell and decay tests, and sine sweep (or 
frequency response) tests. The transfer fimction frequency response data are generally used 
' 20 to identify frequency areas of interest to be investigated fiirther with dwell and decay time 
history methods. The flight time required for the dwell and decay method is extensive 
compared with the flight time required to obtain the frequency response data. Although 
desirable results have been achieved using the previously existing test and analysis methods, 
the complexity of aeroelastic flutter and the continuing advances in aerospace vehicle 

25 requirements are placing increasing demands on those who design and analyze flutter test 
data of aeroelastic structures. Therefore, there is an unmet need for improved methods of 
analyzing flutter test data, particularly those which may provide improved results or which 
may reduce the required amount of flight testing. 

30 Summary of the Invention 

The present invention is directed to methods and systems for analyzing of flutter test 
data using non-linear transfer fimction frequency response fitting. Embodiments of methods 
and systems in accordance with the present invention may advantageously provide improved 
characterization of flutter test data, may provide improved computational efficiency and 
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The avoidance of the unstable condition of flutter and the determination of the 
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reduced costs, and may significantly decrease the amount of flight testing required in 
comparison with prior art methods. 

In one embodiment, analyzing flutter test data includes reading a plurality of data 
points, each data point representing a motion of an aeroelastic structure {e.g. an aircraft 
5 surface) at a different location. A closed form fit to the plurality of data points is performed 
to obtain an initial curve fit condition. At least one non-linear transfer function frequency 
response curve fit is then performed to the plurahty of data points. 

Brief Description of the Drawings 
The preferred and alternative embodiments of the present invention are described in 
10 detail below with reference to the following drawings. 

FIGURE 1 is a flowchart showing a process for fitting a flutter frequency response 
data set using non-linear transfer function frequency response fitting in accordance with an 
embodiment of the present invention; 

FIGURE 2 is a mathematical derivation of a set of analytically-derived sensitivities 
15 for a State Space Model in accordance with an embodiment of the present invention; 

FIGURE 3 is a mathematical derivation of a set of analytically-derived sensitivities 
for a Pole Zero Model in accordance with another embodiment of the present invention; 

FIGURE 4 is a system for analyzing flutter test data in accordance with an 
embodiment of the present invention; and 
20 FIGURE 5 is a correspondence between a representative result of a non-linear 

transfer function curve-fit with a representative set of test data in accordance with an 
embodiment of the present invention. 

Detailed Description of the Invention 
The present invention relates to methods of analyzing of flutter test data using non- 
25 linear transfer function frequency response fitting. Many specific details of certain 
embodiments of the invention are set forth in the following description and in FIGURES 1-5 
to provide a thorough understanding of such embodiments. One skilled in the art, however, 
will understand that the present invention may have additional embodiments, or that the 
present invention may be practiced without several of the details described in the following 
30 description. 

Analyzing flutter test data in accordance with embodiments of the present invention 
may involve non-linear optimization to fit a transfer function frequency response to a set of 
test data to provide an improved characterization of the natural frequency and damping 
characteristics for each mode of the fit. Embodiments in accordance with the present 
35 invention may overcome several disadvantages of the prior art analysis methods while 

■ 3 * BLACK LOWE & GRAHAM^^ 

^ 

t J- BOEI-MI84AP , . ^ 

25315 Second Avenue 

PATWTRADEMAKKOFFICE ^Jf'^^J ^^^^^""1 l^^Lr 

206.381.3300 • F: 206.381.3301 



providing an efficient method to fit the gain and phase of transfer function firequency 
response data. Methods in accordance with the present invention may also allow an analyst 
to increase a number of modes in the fit with relative ease to improve the fit with the 
experimental data, thereby allowing improved fitting capability. 
5 FIGURE 1 is a flowchart showing a process 100 for fitting a flutter time history data 

set using non-linear transfer fimction frequency response fitting in accordance with an 
embodiment of the present invention. In this embodiment, the process 100 begins by reading 
a plurality of test data points at a block 102. In a particular embodiment, the plxirality of test 
data points may represent flutter test data, and each data point may represent an amplitude 

10 value at a given test time. 

At a block 104, the pluraHty of test data points are transformed into a transfer 
fiinction frequency response (TFR). Well established methods are used to compute the TFR 
from the measured test data. As described more fiiUy below, the process 100 shown in 
FIGURE 1 then fits the transformed data with a State Space Model (SSM) as defined by 

15 Equations (1) and (2). A complex response of the SSM is defined by Equation (3), and the 
Gain and Phase of that complex response is computed with Equations (4) and (5), as follows: 



20 



x = Ax + Bu 
y = Cx + Du 
H(s) = C(sI-A)-*B + D 



(1) 
(2) 
(3) 



where 



25 



30 



35 



X = 

u = 

y 

H(s) = 

s = 
I 

B,C,D = 



mode or state matrix 

matrix of system inputs 

matrix of sensor outputs 

transfer fimction frequency response 

Laplace variable 

identity matrix, and 

coefficient matrices 



state matrix of aeroelastic modes 
All 
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Ai,i = aeroelastic modes 



0 1 



0) = frequency of oscillation 

^ = viscous damping 

and where 

G(s) = 201ogio(H(s)) (4) 

10 P(s) = tan"^(Im(H(s)/Re(H(s)) (5) 

where Re(H(s)) = real component of H(s) 

Im(H(s)) = imaginary component of H(s) 

15 With continued reference to FIGURE 1, the process 100 may then perform a linear 

optimization (or closed form) fit of the data at a block 106. In one embodiment, the closed 
form fit is performed using a well-known orthogonal polyreference (OPR) technique. In an 
ahemate embodiment, a known predictor-corrector fit technique may be used. The closed 
form fit of the block 106 provides only an approximate initial fit of the response data. In 

20 some cases, for example, the closed form fit of block 106 may be so approximate that it may 
even predict unstable modes which are known to not exist. 

Using the closed-form fit as an initial condition, a non-linear optimization fit is 
appUed to the State Space Model (SSM) at a block 107. In one embodiment, a fiinction that 
is minimized during the optimization is a sum of the squares of the residuals, where the 

25 residuals are defined by an error in the gain G and phase P fit curves. In one particular aspect 
of the present invention, the non-linear optimization fit of the block 107 may utilize 
analytically-derived sensitivities, as described more fiiUy below. The analytically-derived 
sensitivities may be computed directly from the data, rather than be computed by differences 
using repeated fiinction evaluations. Because the sensitivities are computed directly, the 

30 speed of the convergence of the non-linear optimization may be greatly improved. 

As fiulher shown in FIGURE 1, an evaluation of the fit of the response data is 
performed at a block 108. The evaluation may be performed manually by an analyst viewing 
a comparison plot of the closed form curve fit with the test data on a control monitor, or 
alternately, may be performed in any suitable automated or semi-automated manner. In a 

35 particular embodiment, an error function is computed and evaluated to determine whether the 
curve fit is acceptable and accurate. 
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If the evaluation of the fit is determined to be unsatisfactory at the block 108, a visual 
evaluation of the fit may be performed at a block 110. At this stage, the analyst may 
determine whether to vary any parameters involved in the curve fit routine in an attempt 
improve the curve fit characteristics. A pre-determined number of iterations in the non-linear 
fit are typically set by the analyst, and the convergence of the residuals may be evaluated 
manually. Alternately, an automated convergence criteria could also be used. One a 
converged fit is obtained, a determination of the adequacy of the SSM is made. During the 
determination of block 108, two issues may be addressed. The first issue may be the 
presence of one or more peaks in the resonse data not modeled by the SSM. The other issue 
may be large errors in the Gain (e.g. greater than 4 dB) and/or the Phase (e.g. greater than 8 
degrees). 

Next, an estimate of the firequency and damping of the mode needed for the next non- 
linear fit is performed at a block 112. At a block 114, the analyst may optionally add (or 
remove) aeroelastic modes to (or fi-om) the matrix A in an attempt to improve the fit. When 
a peak is present in the response data and is not modeled by the SSM, the fi-equency for the 
mode is estimated by the fi-equency where the peak appears. A fairly low damping value is 
used for the estimate (e.g, ? = 0.01). When large errors are present, the fi-equency estimate 
for the mode is the fi-equency at wheich the maximum error occurs. In this case, a larger 
damping value is chosen (e.g. ^ = 0.1). The additional mode may be added in such a way as 
to not dramatically change the response of the fit equation, but rather, to incrementally 
expand the fit equation to allow an improved fit with the data during subsequent non-linear 
fit operations (block 107).. 

As fiirther shown in FIGURE 1, after the non-linear optimization fit of the data at the 
block 107, the process 100 may return to the decision block 108 to determine whether the 
curve fit is satisfactory. If the curve fit is not satisfactory, the blocks 110 through 114 and 
1076 may be repeated to attempt to improve the fit. Once a satisfactory fit is achieved, the 
process 100 may proceed to output the results of the curve fit at a block 118, and may then 
continue or terminate as appropriate. 

As mentioned above, the non-linear optimization fit at the block 107 may utiUze 
analytically-derived sensitivities (or partials) of the Gain G and Phase P variables, known as 
the Jacobian matrix A. This matrix may be populated using analytically-derived sensitivities 
that may be computed directly fi:om the data. The Jacobian is a matrix of partial differentials 
of the residuals (row) with respect to each of the design variables (column). The residuals 
may be organized as follows: the gain error (dB) for the first fi-equency then the phase error 
(degrees), repeated for each firequency. The design variables may be organized as follows: 
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the two terms (-o?^ and -^co) from each of the aeroelastic mode blocks (Ai,i) of the state matrix 
A, all terms in the input matrix B, all terms in the output matrix C, and all terms in the feed 
through matrix D. The partials for the PZM non-linear fit are derived much in the same way 
as for the SSM non-linear fits. 
5 For example, FIGURE 2 is a mathematical derivation of a set of analytically-derived 

sensitivities for a State Space Model (SSM) in accordance with an embodiment of the present 
invention. As shown Equations (2-1) and (2-2) in FIGURE 2, the analytically-derived 
sensitivities may be determined by taking the partial derivatives of the Gain G and Phase P 
with respect to the state matrix x, with the Gain G and Phase P from Equations (2) and (3) 
10 above recast in terms of the complex transfer ftinction frequency response Z. Following the 
derivation shown FIGURE 2, the analytically-derived sensitivities for the State Space Model 
(SSM) are provided by Equations (2-17) through (2-30) of FIGURE 2, and may be input into 
the Jacobian matrix to perform the non-linear curve fitting of the test data using the SSM in 
block 107. 

15 An advantage of analyzing flutter test data using analytically-derived sensitivities to 

populate the Jacobian matrix is that computational efficiencies are greatly increased 
compared with altemate methods. For example, populating the Jacobian matrix using 
analytically-derived sensitivities (Equations (2-17) through (2-30)) may significantly reduce 
the computation speed required for each non-linear fit iteration by a factor of 3 compared 

20 with altemate methods of populating the Jacobian matrix which rely upon repeated 
differencing and repeated calls to fimction evaluations. This advantageously improves the 
computational speed and reduces the cost associated with flutter test data analysis. 

In yet another altemate embodiment, the non-linear optimization of the block 107 
may be performed using a Pole Zero Model (PZM) (rather than a SSM). Generally, the PZM 

25 employs a ratio of Laplace polynomials to model a Single Input Single Output (SISO) 
transfer fimction. The polynomials may be expressed as a series of second order blocks 
instead of as a single polynomial. For the PZM method, the residuals are organized as 
follows: the gain error (dB) for the first fi-equency (or mode), then the phase error (degrees), 
repeated for each fi-equency. The design variables are organized as follows: the transfer 

30 fimction gain; the ao term then the ai term for the first pole block, repeated for each pole 
block; the bo term then the bi term for the first zero block, repeated for each zero block. 

FIGURE 3 is a mathematical derivation of a set of analytically-derived sensitivities 
for the Pole Zero Model in accordance with another embodiment of the present invention. 
From FIGURE 3, it will be appreciated that the analytically-derived PZM sensitivities 

35 (Equations (3-6) through (3-15)) for the Jacobian matrix for the PZM non-linear fit are 
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derived much in the same way as for the SSM non-linear fits. Again, the above-noted 
advantages of increased computational efficiency and reduced costs may be achieved. 

FIGURE 4 is a system 400 for analyzing flutter test data in accordance with an 
embodiment of the present invention. Unless otherwise specified below, the components of 
5 the system 400 are of generally-known construction, and will not be described in detail. For 
the sake of brevity, only significant details and aspects of the system 400 will be described. 
As shown in FIGURE 4, in this embodiment, the system 400 includes a computer 402 having 
a central processing unit (CPU) 404 and a memory component 406. The memory component 
406 may include one or more memory modules, such as Random Access Memory (RAM) 

10 modules, Read Only Memory (ROM) modules. Dynamic Random Access Memory (DRAM) 
modules, and any other suitable memory modules. The computer 402 also includes an 
input/oulput (I/O) component 408 that may include a variety of known I/O devices, including 
network connections, video and graphics cards, disk drives or other computer-readable media 
drives, displays, or any other suitable I/O modules. A data bus 410 operatively couples the 

1 5 CPU 404, memory component 406, and the I/O component 408. 

The system 400 embodiment shown in FIGURE 4 further includes a data acquisition 
component 412 operatively coupled to the computer 402. In this embodiment, the data 
acquisition component 412 includes a plurality of transducers 414 that may be arrayed on a 
test article for the acquisition of flutter test data. The data acquisition component 412 is 

20 operatively coupled to the computer 402 via a first conmiunication link 416. 

As fiirther shown in FIGURE 4, the system 400 fiirther includes a control component 
420 having a monitor 422 and a command input device 424 (e.g. a keyboard, an audio-visual 
input device, etc.). A second communication link 418 operatively couples the control 
component 420 to the computer 402. The system 400 also includes an auxiliary output 

25 device 426 coupled to the computer 402 by a third communication link 428. The auxiliary 
output device 426 may include a printer, a compact disk (CD) bumer, a storage device, a 
communication port, or any other desired output device. 

In one aspect, a machine-readable medium may be used to store a set of machine- 
readable instructions (e.g. a computer program) into the computer 402, wherein the machine 

30 -readable instructions embody a method of analyzing flutter test data in accordance with the 
teachings of the present invention. The machine-readable medium may be any type of 
medium which can store data that is readable by the computer 402, including, for example, a 
floppy disk, CD ROM, optical storage disk, magnetic tape, flash memory card, digital video 
disk, RAM, ROM, or any other suitable storage medium. The machine-readable medium, or 

35 the instructions stored thereon, may be temporarily or permanently installed in any desired 
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component of the system 400, including, for example, the I/O component 408, the memory 
component 406, and the auxiliary output device 426. Alternately, the machine-readable 
instructions may be implemented directly into one or more components of the computer 402, 
without the assistance of the machine-readable medium. 
5 In operation, the computer 402 may be configured to perform one or more of the 

aspects of the methods of analyzing flutter test data described above. For example, an 
operator 430 may input a command through the command input device 424 to cause the data 
acquisition component 412 to obtain one or more test data sets read at the block 102 
(FIGURE 1). The test data sets may then be communicated fi'om the data acquisition 

10 component 412 to the computer 402. The computer 402 may be configured to perform the 
methods of analyzing the test data sets. For example, a set of software instructions may be 
stored in the computer 402 (e.g, in the memory component 406) that causes the test data sets 
to be read into the memory component 406 and processed using the CPU 404 in accordance 
with one or more methods of analyzing flutter data, including one or more of the processes 

15 described above with respect to FIGURES 1-3. Altemately, one or more aspects of the 
various processes described above may be implemented in the computer 402 using any 
suitable progranunable or semi-programmable hardware components (e,g. EPROM 
components). 

Results of the analysis in accordance with one or more embodiments of the present 

20 invention may be transmitted via the data bus 410 to the I/O component 408. The results 
may also be transmitted to the control component 420 and to the auxiUary output device 426 
via the second and third commxmications Unks 418 and 428. The operator 430 may view the 
results of the analysis method(s) on the control monitor 422, and may take appropriate action, 
including revising analysis parameters and inputs, and continuing or repeating the one or 

25 more embodiments of analysis methods using different test data as desired. 

FIGURE 5 is a correspondence between a representative result 500 of a non-linear 
transfer fimction curve-fit with a representative set of test data 502 in accordance with an 
embodiment of the present invention. Specifically, FIGURE 5 shows a first graph 510 of 
gain G versus fi-equency, and a second graph 520 of phase P versus firequency, in which a 

30 calculated transfer fimction firequency response (TFR) 512 and the frequency response of a 
calculated SSM 514 are plotted in comparison with the test data 502. Generally, a very good 
fit between the calculated gain G and phase P values and the test data 502 may be noted. 

Methods and systems in accordance with the present invention may provide several 
advantages over the prior art. For example, in the embodiments according to the present 

35 invention, the non-linear transfer fimction curve fit is performed in the gain and phase plane. 
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thereby avoiding the skewing that may be experienced using prior art analysis methods, 
which perform the fit in the complex plan and then convert to the gain and phase plane. 
Another advantage is that the analyst may freely increase the number of modes in the fit 
without re-starting the entire analysis procedure. Unlike prior art analysis methods, the 
5 methods in accordance with the present invention enable the analyst to vary the number of 
modes easily and efficiently. Furthermore, the analyst may easily add more modes to 
improve the fit without fear of degrading the fit already achieved. 

Methods in accordance with the present invention may also allow all of the measured 
flight tests data from multiple transducers to be fit simultaneously. This may advantageously 
10 result in modal characteristics (frequencies and dampings) that best fit all of the measured 
test data. 

Yet another advantage is that fit methods in accordance with embodiments of the 
present invention may be superior to other frequency response methods available in 
commercial off the shelf software. The embodiments of the present invention may provide a 

15 very accurate fit with stable roots. Furthermore, the optimization may converge more 
quickly, especially using embodiments of the present invention that employ a Jacobian 
matrix that is populated with analytically computed sensitivities. 

Still another advantage of the embodiments in accordance with the present invention 
is computational efficiency. Unlike other non-linear optimization methods, the embodiments 

20 in accordance with the present invention may exhibit improved efficiency because the need 
for repeated fiinction evaluations required to get the sensitivities of the gain and phase of the 
response due to changes in the fit equation is reduced or eliminated. 

Yet another advantage of analyzing flutter test data in accordance with embodiments 
of the present invention is that the amount of flight testing may be reduced. Using the 

25 disclosed methods of populating the Jacobian matrix using analytically-derived sensitivities, 
the flight testing required to characterize the aeroelastic modes of a structure may be reduced 
to performing sine sweep tests only. Thus, the need for additional flight testing to conduct 
dwell and decay tests to gather flutter test data to populate the Jacobian matrix A may be 
reduced or eliminated, resulting in significant cost savings and faster analysis schedules 

30 compared with altemate methods. 

While specific embodiments of the invention have been illustrated and described 
herein, as noted above, many changes can be made without departing from the spirit and 
scope of the invention. Accordingly, the scope of the invention is not limited by the 
disclosure of the specific embodiments set forth above. Instead, the invention should be 

35 determined entirely by reference to the claims that follow. 
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